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Abstract 

We consider the class of quantum mechanical master equations de- 
fined on a generic Banach space, arising by projecting weakly perturbed 
one-parameter groups of isometries. We show that the possible semi- 
group approximations are far from unique. However, uniqueness can be 
reestablished through the introduction of a dynamical time averaging 
map. The generator of the resulting Contraction Semigroup is always 
well defined, irrespective of the dimensions of the projected subspace, 
and of the spectral properties of its free dynamics. We show how our 
approach includes and generalizes the preexisting literature. 

1 Introduction 

After the pioneering works in 1974 and 1976 by Davies ([H [2]), a huge 
amount of physical information about the irreversibility and the evolution 
of open quantum mechanical systems has been gained [3]. Alicky [3] showed 
in 1977 that these efforts were deeply connected to the celebrated "Fermi's 
Golden Rule", that now had become mathematically consistent. The con- 
ceptual importance of these works is clearly not only academic, as the need 
for a better understanding of irreversible processes has never been more ur- 
gent. Today, so many nanotechnologies are pushing devices towards limits 
where neither quantum phase coherence, nor dissipation/dephasing, can be 
neglected [5j El \7\ EJ [9] . Many attempts to improve the theory have been 
made since then (see for example [10J HI]), but despite the compelling need, 
no substantial, fundamental progress, directly applicable to nowadays tech- 
nologies, has been made so far. 
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To be more specific, the problem is to understand the dynamics of a 
subsystem of interest, when the global system is fully coherent. In many 
cases, the dynamics of the global system can be splitted into a part that 
leaves the subsystem invariant, plus an interaction between the subsystem 
and the remaining "unobserved" degrees of freedom. The problem then 
arises whether or not the subsystem can be given a markovian, possibly 
dissipative, dynamics as a consequence of its interaction with those degrees 
of freedom. It has been shown since the '70s that this problem can be 
attacked, and partially solved [HE], when the amount of unobserved degrees 
of freedom is huge, and the interaction is made small. This last condition is 
referred to as the "weak-coupling limit". 

In [1] the author was able to give a solid physical model of a discrete 
"atom" (system A) interacting with a fermionic particle reservoir at ther- 
mal equilibrium (environment B). In that case, the subsystem was made of 
unentangled pairs of atomic states (also referred to as "density matrices") 
and a fermionic thermal equilibrium state, while all the entangled pairs just 
constituted the remaining, uninteresting, degrees of freedom. The model 
was of high conceptual importance, as the markovian (and dissipative) dy- 
namics for the subsystem was shown to guarantee the state positivity at all 
times. This fact gave just enough internal robustness to the model as to be 
of invaluable practical use, for at any time, the subsystem evolution could be 
given a strong physical meaning. But the system A had to be finite dimen- 
sional, or at least, its unperturbed hamiltonian was forced to have discrete 
spectrum. This fact constitutes a severe limitation with respect to nowadays 
needs to explore systems at mesoscopic scale, where energy spectra are often 
of mixed nature. 

Unfortunately, positivity through the weak-coupling limit procedure be- 
comes no more available when the system is infinite dimensional, or the 
spectrum is continuous. In that case indeed Davies showed that a marko- 
vian approximation, for the subsystem exact projected dynamics, could still 
be achieved [2], but positivity could not be shown even for the old and "safe" 
partial tracing over the thermal bath's degrees of freedom. Being more gen- 
eral, the theory was and is currently applied in many physical urgent contests 
(see for example |12j). but certainly it did not share the enormous success of 
the previous one among physicists, precisely because of the serious positivity 
limitation. For example, all of the steady state analysis became, physically, 
completely meaningless. 

In this work we provide a rigorous solution to this problem in the general 
contest of Banach spaces, by extending and putting on more firm grounds 
our initial study in [13]. We suppose the global system is undergoing a 
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fully coherent evolution (according to a perturbed one-parameter group of 
isometries) , and discover a whole class of new markovian approximations 
for the subsystem, by suitably manipulating the memory kernel of the exact 
projected evolution, in the limit of small perturbation. This class includes 
the markovian approximation performed by Davies in [2J, and every gener- 
ator in it is well defined irrespective of the subsystem spectral properties 
and dimensions. Since we will work in the general contest of Banach spaces, 
we will not address positivity at this time, but it will be evident that none 
of generators found so far gives rise to a contraction semigroup, which is 
a basic requirement for a Quantum Dynamical Semigroup [2] in the more 
specific contest of W*-algebras. 

At this point, we will introduce a "transition time" function, that will 
scale with the inverse of the perturbation, and will represent physically the 
transition time among the set of relevant states, due to the perturbation. 
This time will serve us to perform a " dynamical average" among all the gen- 
erators previously found. The result will be a "dynamical time averaging 
map", very similar to the time average proposed in pQ, but fundamentally 
different in that i) it will scale with the coupling constant and ii) it averages 
among different generators (in pQ the time average acts on one generator 
only). We shall be able to prove that the resulting transition time dy- 
namically averaged generator i) correctly approximates the exact projected 
dynamics in the weak-coupling limit, ii) is always well defined, irrespective 
of the subsystem spectral properties and dimensions, Hi) accounts also for 
first order contributions, iv) boils down to the averaged generator in pQ in 
case of discrete spectrum, and we shall prove the all important result that v) 
it gives rise to a contraction semigroup on the projected Banach subspace. 

Thus the requirement of a contraction semigroup essentially removes the 
degeneracy of all the possible semigroup approximations. The generality 
with which we shall be able to make our statements (no dependence on 
spectral properties or dimensions of the projected subspace) establishes on 
firm grounds the possibility of many new important physical (and mathe- 
matical) applications. 

2 General Theory and Motivation 

We report from Davies [2] the general framework we'll be involved with. We 
suppose that Pq is a linear projection on a Banach space B (that represents 
some global system), put P\ = 1 — Pq and Bi = P{B, so that 

B = B Q ®B 1 , (1) 
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and we take Bq to represent a subsystem of interest, B\ being the remaining 
degrees of freedom. We suppose that Z is the (densely defined) generator 
of a strongly continuous one-parameter group of isometries Ut on B with 

U t P = P U t (2) 

for all t S R, or equivalently 

[Z,P }=0 (3) 

and put Zi = P- L Z. We suppose that A is a bounded perturbation of Z 
and put Aij = P{ APj. We let U x be the one parameter group generated by 
(Z + XA 00 + XA U ) so that 

U t x P = P U t x (4) 

for all t G R, and let V t x be the one parameter group generated by Z + XA, 
so that 



V t x = U t x + X f U t x _ s (A i + A 10 )V x ds. 
Jo 



(5) 



Then putting 

X x = P U X (6) 
and defining the projected evolution as 

W x = P V X P (7) 

and one obtains the all important closed and exact integral equation 

W x = X x + A 2 f ds I" du Xl s A 01 U x _ u A w W x . (8) 
Jo Jo 

This is nothing but the integrated form of the well known master equation 
constructed by Nakajima, Prigogine, Resibois, and Zwanzig [151 116] , which 
follows by repeatedly making use of ([5]) with 

W x = X X + X f ds X X _ S A 01 V X P (9) 
Jo 

and 

PiV/Po = A f du U X _ U A 10 W X . (10) 
Jo 

Now let X t A is a one parameter group of isometries. For example, Davies 
proves [2] that 
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Lemma 3 If \\Pq\\ = 1 and ||e || = 1 for all t £ I, then X* is a one 
parameter group of isometries on Bq for all real A. 

Then, changing variables to x = s — u, a = X 2 u and introducing the 
time rescaled (and j4oo- r enormalized) interaction picture evolution Wr' 1 = 
X_ X -2 T W^2 T , one is led to 0: 

= 1 + f da X* x . 2a K(\, r - a)X^ 2a Wy (11) 
J o 

where 

/•A" 2 T 

K(X,t)= / X_ X A 01 U£A 10 dx. (12) 
J o 

This form separates an "interacting" and "slowly varying" part K(\, r) from 
the "rapidly oscillating" free-evolution X x -2 a to second order in the coupling 
constant A. Now in the weak-coupling limit A — > 0, the slowly varying 
integral kernel K(A, r) converges to 

/•oo 

K D = / U^ X A 01 U X A 10 dx , (13) 
Jo 

where Kd is the celebrated Davies' superoperator. Substituting Krj ~ 
K(X,t) in (fTT|) and moving back to the "Schrodinger picture" we obtain 
the markovian approximation for our subsystem dynamics 

Indeed in [2j an important theorem shows that under reasonable and general 
conditions the approximation holds in the weak-coupling limit, up to A~ 2 - 
rescaled positive times r, independently of the subsystem dimensions or 
spectral properties. Unfortunately, Kd does not guarantee positivity of the 
generated evolution, for example in the case of partial tracing over a bath 
PQ, and thus may fail to furnish a physically acceptable model at times 
t > A~ 2 r. This fact precludes the important possibility to study the limit 
dynamics at all times, together with all the steady state analysis. 

However, in [1] the author shows that, in case Bq is finite dimensional (or 
more generally when Zq has discrete spectrum), one can define a temporal 
average K^ D of and show that the generated evolution is still asymptotic 
to the exact one in the weak-limit (at least when Aqq = 0). Moreover, the 
author also shows that if the projection Pq is taken to be the partial trace 
over a bath, the averaged K D generates a Quantum Dynamical Semigroup 
(QDS) in the sense of p[4j. 
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So finally Kjy does not guarantee a positive evolution at all times, while 
Kp is well denned only in the far too restrictive case of discrete Zo-spectrum. 
This serious problem is physically so important as to motivate our work. 



4 Classes of Generators and Dynamical Transition 
Time 

In our first main theorem, we shall look for the most general form of an 
operator whose evolution is still asymptotic to the exact one in the weak- 
limit. To this purpose, note that Kq has been defined thanks to the change 
of variable 

J" A" 2 [ x 



-2 n I I ~ I ■ ( 15 ) 



Here instead we would like to allow for the most general change of variable 
that keeps a A~ 2 jacobian, proper of a second order approximation, while 
fixing the relative variable to be s — u = x. So we chose 



s \ _ ( A 2 a + \ \ ( a - X 2 q 



u ) V A 2 a — \ ) \ x 



(16) 



with inverse 
a 



., v (l "" )A! (i+ i A2 )(:Mt). 

for some real a and q. Some straightforward algebra shows that the in- 
tegration domain s = . . . A _2 r, u = 0...s in ([8j), becomes the domain 
P(A, r, a, q) in the (a, x) plane given by the triangle of vertices 

V(X, r, a, q) = A | (X 2 q, 0), (r + X 2 q, 0) Y Q - a) r + X 2 q, X' 2 ^ X (18) 

(see figure (P)). Accordingly, equation ([8]), written for W T = X_^-2 T W£- a , 
becomes 

W*'* = 1+ I! dodxX x ,„ A 01 U*A 10 X* , n 



T>(\,T,a,q) 

xW*'\ 2{q+{a „ 1/2)x) . (19) 



We now consider the following facts, which will be made precise and proven 
in Theorem l7.lt 
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X4- 




ftJSfll S>(\T,a,q) 



+X a q (a-J4)T+A a q 0" 



Figure 1: Integration domains T>(X,T,a,q) and D(A,r, a, q) denned in (fT8|) 
and (f32l) respectively (we have put q > , a < — 1/2 for clarity). The arrow 
indicates the asymptotic behavior of the two domains in the weak coupling 
limit A —* (see discussion in the text following (|18() ). 



T>(\, t, a, q) — ► [0, r] x [0, oo) as A — > 0, for any real a and g. This 
justifies the approximation 



da dx 

V(X,r,a,q) 



da dxe-^ 2 / T ^ 2 

J 



A ss 0, (20) 



for some real function T of the coupling constant, provided lim^o T(X) 
+oo. Indeed, exponential factor would not change the asymptotic be- 
havior of the kernel, and would naturally conserve the original char- 
acter of a bounded integration domain for every nonzero A; 



In the weak coupling limit A — > one could approximate 



W 



X.i 



cr-X 2 (q+(a-l/2)x) 

in the integral kernel of (119p . 
Once all this is proven to be legitimate, equation (|19f) becomes 



(21) 



XA 



1 + 



da X 



x 

-x- 



(a,q,T(X))X X -2 a 



X 



w. 



X.i 



or, which is the same, 



W t x « exp{(Z + A^oo + A 2 if (cM?iT(A)) )t}, 



(22) 
(23) 
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where we give the following 

Definition 5 Let a, q € R be two given real numbers, and let T G C{[— 1, 1],M 
be a real valued positive continuous function on the interval [—1,1], where 
M = lU {°o}. For X £ I define the linear operator Kr a> qT(x\) on &o as 

K (a , q , n x)) = / dx e -(f ) / T W 2 Wfl A)l^iotff a _jv_ fl . (24) 



Denote also with 

W\ = exp{ (Z + AA 00 + A 2 ^ (Qj(?iT(a)) )t} (25) 
i/ie associated semigroup on Bo- 

Then, under fairly general hypotheses, we shall prove that is compatible 
with the exact in the weak coupling limit A — > 0, up to A _2 -rescaled 
positive times r: 

lim i sup ||W t A - W?|| I = 0. (26) 
[o<t<\- 2 - J 

We stress that, as in [2], we do not assume that B$ is finite-dimensional. 

The following lemma is not new, and for example it is contained in 
Theorem 1.2 of [2], but we report it here as we shall make use of it repeatedly 
all throughout. 

Lemma 6 Let b G Bo be given, together with some real f > 0. Suppose 
and w\ are operators on Bo such that f\(r) = X\_ 2 Wv_ a 6 satisfies 



n>0 

for a Volterra operator Tlx on the Banach space V = C°([0,t],B ) of contin- 
uous Bo-valued functions on the interval [0,r] (assume the same holds also 
for W t , with associated / A and Ti\). Suppose there exists a real positive c 
such that ||H A || V < cr and \\Ti.\\\ v < cr uniformly on |A| < 1. Put 

i) liniA^o \\H\ -H\\\v = ' 

it) lim A _> ||/ A - /J,*, = ; (28) 
Hi) lim A ^ sup < t < A -2 T ||W t A 6 - W t 6|| Bq = . 

Then i) ii) Hi). 
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Proof. Of course ii) =>- in) as is group of isometries and one has 

sup ||W^6-W?6|| Bo = sup ||/a(t)-7aWIIb = II/a-7aIIoo. (29) 

0<t<\- 2 T 0<T<T 

Then, subtracting the von Newmann expansions for f x and f x one obtains 

oo 

||/a-7aIIoo<^II^-^6|| 00 

71=1 

OO 

= E H^a 6 - n V^ h + ^T 17R xb -■■■ + n x nT l b - W^Hoo 

n=l 

oo 

= E WKrHnx - n x )b + n r x l - 2 (n x - n x )n x b + ... + (n x - Ha)^ _1 &I|oo 

71=1 

< \\nx-n x \\jb\\ Bo J2\^ (30) 

?1=1 ^ ' 

and the last series is (obviously) convergent and independent of A. Note 
that we have used (and will use throughout) the important property that if 
H is Volterra and \\H\\ < C, then \\H\\ n < C n /nl. 

This shows that i) => ii) and thus finishes the proof. □ 

We still need a technical but important result, that will allow us to per- 
form approximation (|2 1 [) : its interpretation will become clear in the context 
of Theorem 17. 1\ but it deserves to be reported in the more autonomous en- 
vironment of a Lemma, as it will find application also in our second main 
result, Theorem 19.11 

Lemma 7 Let Bo be a Banach space, T > 0, and let V = C°([0,t],Bq) be 
the Banach space of continuous functions from [0,r] into Bq. For some real 
q and a, letT>(X,T,a,q) be the triangle in the (o~,x)-plane of vertices 

V(X, r, a, q) = A ( (X 2 q, 0), (r + X 2 q, 0), ( ( ~ - a ) r + X 2 q, X^r] ) (31) 



, , 2 

and define the truncated domain ( see figure (fTJ) ) 

V(X, r, a, q) = V(X, r, a, q) n [0, r] x [0, oo). (32) 

Let T G C([— 1,1], R) be a real valued positive continuous function on the 
interval [—1,1], and assume that T(X) ~ |A| ? for X -» with £ < 2. Let 
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7~i(\aq) be a Volterra integral operator on V and assume it can be put in the 
form 

(n {Xaq )9){r) = jj dadx e" W K aq (\, a, x) g (a - A 2 (q + (a - 1/2) x)) , 

V{X,r,a,q) 

(33) 

for a suitable kernel K aq (X,a,x), and define the Volterra integral operator 



(n {Xaq) g)(T) = JJ dadxe W K aq (X, a, x) g{a). (34) 

V{\,T,a,q) 

Suppose that \\K aq (\, a, x)\\ = k(x) independently on a,X,a,q, and that 

i 

dx k(x) = c (35) 

for some finite < c < oo, and assume WHfXaq)\\ < 

Let b £ Bo and define the von Neumann series f(x a q) an d f(\aq) through 
equation [27\ ) in Lemma 

Then __ 

WH(\ aq )\\ < tc (36) 

and 

lim Wf(Xaq) ~ l(Xaq) Woo = . (37) 

Proof. Estimation (|36|) is a trivial consequence of the definition of 



\\n {Xaq )9(r)\\<T fcO*ONU. (38) 
Jo 

To show the validity of (I37p we start proceeding as in Lemma [6] to obtain 

ll/(Aag) - f(Xaq)\\oo < /2 z2 („ _ Al ^{Xaq) ~ H(Aag)) ^{Xaq) b \\oo ■ 
n=l 2=1 { h 

Note that the case I = has been dropped, since H(\ aq )b = H(\ aq )b trivially, 
as b is constant, as one can see from (]33[) . Now if we could show that for 
every e > there exists some A > such that |A| < A implies that for every 
1>1 

— —l (re) 1 - 1 _ 

\\(7~L(\ aq ) ~ ri(\ aq) ) W(Aa 9 ) &||oo < ^ _ ^, r ||&||b e , (39) 
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we would be done, as, following estimation (f39]) . we would have 



n-1 



Wh-f 



Alloo 



< 



EE 

71=1 1 = 1 



(tc)"- 2 
(n-l)\(l-l)\ 



I So 



(40) 



and the series would obviously converge (because c > by hypothesis). To 




Figure 2: Integration domains V(X,a,a,q) (rectangular triangle (1)), 
T>(\, a — X 2 (q + (a — |) x) ,a,q) (rectangular triangle (2)) and S(X, a, x) 
(region (1) U (2)). We have put a < —1/2, q > for clarity. The two 
truncated triangles T> are similar, as are also the untruncated, underlying 
T>. Note the scaling behavior with A 2 . 
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show that property (|39|) holds, we take I > and evaluate 



If 



dadxe T < A ) 2 K aq (X,a, x 



x {(n\ Xaq) b) {a - X 2 (q + (a- l/2)x)) - {n\ Xaq) b){a)) 



(*/2) Z 



dadx e T < A ' 2 K aq (X,a,x) 



V(\,r,a,q) 



daidxx A A (ai,a;i)e K aq (X, a h Xl ) (H {Xaq) b)(a$l) 

S(\,a,q,cr,x) 

where we have put A A = Xi ~ x\i Xi being the characteristic function of 
V(X, a—X 2 (q + (a — |) x) ,a,q), X2 the characteristic function of V(X, a, a, q), 
and we have defined (see Figure [2|) 



S(X, a, q, a, x) = V ( A, a — X 2 ( q + ( a — - ) x ) , a, q I U V(X, <7, a, q) . 



(42) 



Passing to the norms, we estimate the left hand side of (139H as 



W2) Z 



dadxe T W 2 k{x) 



< sup 

0<T<¥_ 

V(\,T,a,q) S(\,a,q,a,x) 



da\dx\ |A (o"x, x\)\ k{x\) 



< 



(tc) 1 - 1 



(*/2)^ 



..II&Hbo sup [[dadxe k(x) E(X,a,x) 

— -LJ! o<r<- J; 



©(A,T,a,g) 



0-1)! 



0<o-<r JO 



< v . ; ... \\b\\ Bo T sup / die n^k(x)E aq (X,a,x) 



1~L(\ aq) b 



(43) 



where we named 



i aq (X,a,x)= jj da 1 dx 1 \A x (a 1 ,x 1 )\k(x 1 ) 



S(\,a,q,(T,x) 



(44) 
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At this point let us consider the sector defined by q > and a < —1/2. In 
this case, it is but a straightforward algebra to show that the two triangles 
T> (A, a — X 2 [q + (a — ^) x) , a, q) and V(X, a, a, q) (the same applies to V) 
are similar, and that their left edges lie on the same line of equation 

x l (X,a 1 ) = 1 (a x - X 2 q) (45) 



a 



in the (o"i, x\) -plane (see figure ([2])). In particular, 



For this reason we easily compute 

1 



2?(A,(7,a,( 



> 

< l 



E aq (X,a,x) 



< A 2 



< cA z | 9 + 



sign 



1 



a x 

2 



a-\ 2 (q+(a-±)x) /•a: ! (A,<7 1 ) 



do - ! 



+ I a - - | x 



1 

°"2 



d^i k(x\) 



dxi k(xi] 



(46) 



According to the last line in estimation (|43j) . we study 



sup 

0<<r<r Jo 



die T ( A ) 2 /c(x) H ag (A,o-,x) 



< c 2 A 2 g + c 



1 




a 




2 





(z/2) 2 



dxk(x)e iwA 2 x 



(47) 



The first term goes to zero with velocity ~ X 2 q, whereas the dominated 
convergence theorem applies to the second term, showing convergence to 
zero uniformly on q, with velocity ~ A 2 ^^ (note that we have supposed 
£ < 2 strictly). This, when put in (f4*3j) . shows the validity of (|39|) . 



A very similar analysis can be done for the remaining sectors q ^ 0, 



> 



\a\ < 1/2 and a > 1/2. In each case, the net result is that the estimation 
(|4T|) is always of the form 



sup 

0<cr<T Jo 



(*/2) 



dxe ^^(x) E aq (X,a,x) < dX 2 \q\ + C 2 X 2 ^ 



(48) 



for suitable real positive constants C\ and Again, this shows the validity 
of ([39]) . and thus concludes the proof. □ 



We are now in position to state our first main result: 
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Theorem 7.1 Suppose that is a one-parameter group of isometries. 

Let T G C([— 1, 1], R) be a real valued positive continuous function on the 
interval [—1,1], and assume 

T(A) ~ |A| - *f, A^O (49) 

for some real positive T > and < £ < 2 (strictly). 

Suppose that there exists some < c < oo suc/i i/tcrf /or every r > 

/ Poi^^Mioll^ < c (50) 

JO 

is bounded uniformly on |A| < 1. 

Suppose also that for every < r < oo 

lim / 11^01(^-^)^1011^ = 0. (51) 
T/ien /or every T > 



lim { sup ||W t A - W7|| } = 0. (52) 
[o<i<A- 2 r J 

Proof. Let V be the Banach space of norm continuous i3o-valued func- 
tions on [0,r], and let b £ Bo. Define the "interaction picture" time rescaled 
solution of © 

f x (r) = X\^ T W^ 2r b. (53) 
Then fx is a solution to the integral equation 

fx = b + H x fx, (54) 

where the integral operator Tlx is defined (recall that is a group of 
isometries) by 

(Hxg)(r) = A 2 [ X ds f du X*A 01 U*_ U A W X* g{\ 2 u) . (55) 
Jo Jo 

In [2] Davies shows that 7i\ is a Volterra operator. Indeed, by changing 
coordinates according to (fl~5]) . Eq. ([55]) can be given an explicit Volterra 
form, namely as 

(Hx 9 )(t) = [ T daX\_ 2a K(\,r- a)X^ 2a g(a) , (56) 
Jo 
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where we defined the "slowly varying" kernel 

r \- 2 r 

K(X,t)= dx X* x A xU£A w . (57) 

Jo 

Because of this reason, and since K(\, t) is manifestly bounded by c, uni- 
formly on A (thanks our boundedness hypothesis), it follows that 

\m\\ < c n r n /n\, (58) 

and also that the associated von Newmann series expansion 

f x = b + n x b + n 2 x b + --- (59) 

converges. 

We can proceed in similar fashion also for the semigroup (|25h : iteration 
gives indeed 

WUr = X\-*r + [ ^ X^. 2iT _ a) K (a!q>T{x)) W^ 2a . (60) 
Accordingly, we define 

J x (r)=X\^ T WUrb (61) 
so that it follows that f x is a solution to the integral equation 

7x = b + n { x aq) 7x, (62) 

where we have defined 

(n iXaq) g)(r)= [ da X^ x . 2a K {a ^ T{x)) X^ 2(T g(a). (63) 
Jo 

Now again, Ti.( Xaq ) is a Volterra operator, and since < c, equations 

([58]) and ([59]) follow analogously for TL( Xaq ). 

Clearly (see Lemma [6]), we must show that for any chosen b £ Bq, 

sup \\W t x b-W}b\\ = ||/a-7aIU -0, A^O. (64) 

0<t<A" 2 r 



We shall do that by defining suitable Volterra operators, denoted with 
n< (\aq)i 3 = 0...N, such that nf x \ q) = H\, WjJ^ = H {Xaq) , and either 
we can show, according to Lemma [6l that 

iW&j-WfcSllv-O, A-0 (65) 
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or more directly that 

ll/( ( £ 9 )-/(t g )IU-0, A^O (66) 

where 

oo 

/ ( t) = E« a ,)"' m 

n=0 

is the associated von Neumann series. Then, our conclusion (see Lemma [6|) 
would follow fro 

N 

II/a-JaII^EII/m-zSii- (68) 

J'=l 

To follow our purpose, instead of ()15p . we perform the coordinate trans- 
formations (|16p As we noted before, the integration domain s = . . . A _2 r, 
u = . . . s, becomes the domain T>(\, r, a, g) in the (a, x)-plane defined in 
(fT8|) . and depicted in figure ([T]). Accordingly, the integral kernel in ([55]) is 
now written as 

(^)(r) = // ^^ A - 2CT+9 -(a + l)^oi^io4_ v _ ?+(a _ |):c 

V(\,T,ct,q) 

xg(a -X 2 (q + (a - l/2)x)) . (69) 

We shall work throughout with the choice a < —1/2 and q > 0, without 
loss of generality (the structure of the proof is basically the same for all 
the remaining sectors, which we leave to the reader, and actually the case 
\a\ < 1/2 presents less difficulties). 

The first thing we shall be concerned with, is to find a way to substitute 
the free "polarization evolution" P\U X in place of the interacting P\U X in 
the middle of the kernel in (|69p . This step is accomplished by defining a 
related integral operator 7iy^ by 

^(x aq )9)(r) = /J dadx X\_ 2u+q _ { ^ )x A Q1 U x A 10 Xl_^_ q+{a _ l)x 

25(1) i\T,a,q) 

x g( a -\ 2 {q+ {a -l/2)x)) . (70) 

Note that we have attached the subscript "(Aag)" all throughout: although some 
Volterra operator may not actually depend on a, nor q, this unifying notation will become 
useful in the sequel. 
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Here we have denoted T>i 



T> for sake of notation. Note that TC 



(i) 



does 



'(1) - - *™ ^ ^« uc "-(Aog) 

not depend on a, nor on q, although both the latter parameters appear in 

its definition. To compare the two integral operators, take a bounded g £ V 

and estimate 



TL\ — TL 



(i) 

(Xaq) 



9 



< 



dadx \\A i(U* - U X )A W \\ \\g\ 



T>(\,T,a,q) 

,A" 2 r 



< max{l, |l/2 — a\} r\\g\ 



dx \\A 01 (U£-U X )A 



10 I 



(71) 



o 



(note for later purposes that the estimation does not depend on q). Our 
convergence hypothesis (|5ip on U x allows then to conclude that this goes to 
zero when A — > uniformly on every \\g\\ = 1, so that we obtain 



lim 

A^O 



TL\ — TL 



(i) 

(Xaq) 







(72) 



We proceed along similar lines to smooth the kernel with T(A): define 



dadx e W X 



A 



A 



A -2 CT+g _( a+ i)^01^10^ A - 2(J _ g+ ( a _ 

f (2) (A,T,a!,g) 

x 5 (<7-A 2 (<z + (a-l/2)x)). 

Here we have denoted V( 2 ) = ^(i)- To compare with %, take a bounded 
<7 G V and estimate 



H 



(i) 

(Aag) 



(2) \ 
(Xaq)) 



< 



dadx\l-e-^ 2 / T ^ 2 \ 



\A m U x Ai \\ \\g\\ 







1 


1 


< max < 


a 








2 





T>(X,T,a,q) 



(^/2) 2 



dx|l-e \\A 01 U X A 10 \ 



o 



(74) 



Hypothesis ([50]) furnishes an integrable upper bound to the last integrand, 
and so the integral goes to zero in the limit A — > because of the dominated 
convergence theorem: in fact, one has pointwise convergence 



1 - e T W 



0, A 







(75) 



due to our hypothesis < £. Uniform convergence on all \\g\\ = 1 in the last 
line of our estimation then shows that 



lim 

A^o 



n 



(1) 

(Xaq) 



u 



(2) 

(Xaq) 







(76) 



17 



(2) 

(Xaq) 1 



as T^xaq) an d T~t\, is also a Volterra operator, as is 



We note that H 

easy to verify: it would suffice to apply the inverse transform (|17p for the 
specific choice of a and q, and subsequently apply the transform (fT6j) for 
q = and a 

find that it does not depend on (ct,q)) 

Now define the restricted domain (see figure dU) 



(2) 

1/2, to find 'HWaq) ^ s explicit Volterra form (and also 



V {3) (X,T,a,q) =V {2) (X,T,a,q) n [0,r]x[0,oo) 



(77) 



(3) (2) 

and define Ti^ aq ) accordingly, as the restriction of Ti)^ a s to V^(X, r, a, 
that is, 



_MH)1 
dadx e X x 



-\~ 2 a+q-(a+\)x 

v (S){\r,a,q) 

x g( a -\ 2 {q + {a-l/2)x)) . 



i\AoiU x AioX 



A 2 <r—q+(a 



Of course this is again a Volterra operator, as the image of T>r 3 \(X,T,a,q) 
under the composition of (fT7|) and (fT6|) . for g = and a = 1/2, is inside the 
rectangular triangle, image of T>(2)(X,T,a,q) through the same transforma- 
tions. So the composition of (fT7|) and (fT6j) . for q = and a = 1/2, will put 



both n 



(2) 



(3) 

(Aag) anc ^ ^(Xaq) * n their explicit Volterra form. 
We must prove that 



lim 



Ti 



(2) 

(Xaq) 



Ti 



(3) 

(Xaq) 







(79) 



/(2) 



(note that for \a\ < 1/2 this is obvious, as one has 7~t( Xaq ) 
upper vertex of the triangular domain projects on the triangle cr-axis basis, 
and thus £W) (A, r, a, q) = T>r 3 \(X,r,a,q)). To this end, we take as usual 
any g £ V with \\g\\ = 1 and estimate 



Here 



(2) -,(3) 



n\ ' - n 



< h\\ 



(Xaq) 



9)(r) 



I I dadx II-AoiE^AloII \\g\\oo 

V(X,T,a,q)n{r<a<\l/2~a\T+X 2 q} 



poo p(^-a)T+X 2 q 

1 dx IIAji^Aioll + 



o 



da / dx \\A iU x A w 

X 2 q Jx r (X,q,a) 



x r (X, a) 



q + X-\r-a) 
a + i 
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is the rr-coordinate of the right edge of the triangle T>(\, r, a, q) (see figure 
(H|)). The first term in the curly brackets clearly goes to zero uniformly on 
r (with speed ~ A _2 g), because of our boundedness hypothesis (f50|) . In 
the second term in the curly brackets, we change coordinates according to 
(j — ► (J — X 2 q, so that it becomes equal to 

/ da dx \\A 01 U X A W \\. (82) 

Jt J x r (\,0,(j) 

Now since x r (A, 0, a) > a.e. in the a- integration domain (everywhere 
except a = r), we have the following convergence (pointwise with respect to 
o) 

r+oo 

lim / dx \\A 01 U X A W \\ = , (83) 

thanks to hypothesis (|50p for A = 0. Since the integration domain for a 
is bounded, convergence to zero follows for the first integral in the curly 
brackets in the estimation above, because of the dominated convergence 
theorem. Uniform convergence on all < r < r (and real q) follows from 
the fact that the u-integration domain is compact. This proves convergence 

According to our roadmap, we now define the following "time localized" 
Volterra operator: 

i^ aq) 9)(r) =JJ dadxe^ X^^^^A^^^ 

(84) 

where we have put P( 4 ) (A, r, a,q)) = (A, r, a, q). It turns out that prov- 

(3) 

ing an operator convergence to j is impossible, due to the strong re- 
quirement of uniform convergence with respect to any normalized g G V. 
However, the Volterra operators 'H^ aq \ and 'H^aq) ^° ^ m ^ u * ne hypotheses 
of Lemma [71 and so we conclude that 

We now go back to consider the domain T>u\ (A, r, a, q) as a function of 
A, and note that it tends to fill the strip 

Z>( 5 )(t) = [0,t]x[0,oo]. (86) 
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Accordingly, we define the following Volterra integral operator on this strip: 



H {5) a 
v ' Jo Jo 

and note that 



da I dxe X x 



AmUtAmX 



x 

' \~ 2 o—q-\-(a- 



(87) 



n 



(5) 

(\aq) 



n 



(4) 

(Xaq) 



9)(r) 



< 



j j dadx \\AqiU x Axq\\ \\g\ 



D(5) (r)\n( 4 ) (X,T,a,q) 



< 



where 



dx ||A)if^4io|| + 



xi(X,q,a) 



da 



X 2 q 



\~ 2 a 



dx \\A iU x Aii 



xi(X,q,a) 



(89) 



is the x-coordinate of the left edge of the triangular domain V(\, r, a, q) (see 
figure dl])). Now the first term in the curly brackets clearly goes to zero as 
A — ► (with velocity ~ A 2 g), due to hypothesis (f50|) for A = 0. In the second 
one, as before, we change coordinate according to a — > a — X 2 q, obtaining 



[T-X 2 q r 

/ da 

JO Jx 



dx \\A 01 U X A W \\ < 



xi(X,0,a) 



da 



dx \\A 01 U X A W \ 



(90) 



This last term can be seen to go to zero uniformly on < r < f (and real q) 
by invoking, as done before, the dominated convergence theorem, and using 
the fact that x~i(\,0,a) > a.e. in the (compact) cr-integration domain 
(everywhere except a = 0). So it follows that 



n 



(4) 

(Xaq) 



n 



Now we can finally compare with Ji 
ing obvious terms, we estimate 



(5) 

(Xaq) 
(6) 



(Xaq) 



0, A -» 0. (91) 
Tt(\ aq )'- by adding and subtract- 



Ti 



{ (Xaq) - W (Aa 9 )) 9 
oo 

dx \\X X / in -P U 

oo 

X 



\A 01 U x Ai \\\\g\ 



o 



o 



dx ll^oiC/^ollllX^^-PoC/^l^llll.'/lix 



(92) 
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Again, uniform convergence to zero follows by hypothesis (|50p together with 
the dominated convergence theorem, using the fact that for every x £ R 

lim \\X^-P U X \\ =0. (93) 

A^O 

This proves the estimation in (|68p and thus finishes the proof. □ 



7.1 Comments to the theorem 

• The results of the theorem clearly generalize [2] , as the semigroup stud- 
ied there is generated by our particular choice Kd = -f^( Q= i/2, g =o,T(A)=+oo)- 
It would in fact be possible to prove the validity of Theorem 17.11 
also for the case T(A) = +oo identically, if the particular choices 

a = ±1/2, q = are made. 

• For each choice of a, q and T, the corresponding generator Zq + \Aqq + 
K(a,q,T(\)) is always well defined for all A G [—1,1], A / 0, no matter 
which are Bq dimensions or Zq's spectral properties. 

• The function T : [—1, 1] — > R + can be considered as a measure of 
the system's transition times. Indeed, an obvious and natural choice 
would be to put 

r(A) = WW (94) 

With this choice, it is interesting to note that if we chose r = ||^4|| _1 
in the thesis of Theorem 17. H the latter can be rewritten as 

lim J sup \\W t x - Wj || > = 0, (95) 

[o<t<A-!T(A) J 

so that we see that the approximation is valid up to the "dynami- 
cal observation time" A _1 T(A), which is greater than the "transition 
time" T(A), but shorter than possible Poincare-like recurrence times, 
as A m1 T(A) < +oo for any nonzero A. This scaling transition time 
T(A) will play an important role in our second main theorem, through 
the definition of a dynamical time average. 

• Even if the generator in Theorem l7.1l is always well defined for nonzero 
values of the coupling constant A, it does not give rise to a Contraction 
Semigroup, as can be easily seen. This means that we still have to 
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perform some kind of temporal average, like the one introduced in pQ, 
in order to obtain a contractive semigroup. We shall do it in the next 
section. 



7.2 A Sufficient Condition For the Hypotheses 

We provide a sufficient condition for the validity of the hypotheses in Theo- 
rem [77TJ which is a slight adaptation of Theorem 1.3 in [2J, and is physically 
supported by perturbation argumentations. First, define 

ft rtn-l 

an(t)= dto--- dt n \\AoiU to - tl A u U tl -t 2 A u . . . AnU^AwW, (96) 
Jo Jo 

which are, as one can easily realize, coming from the expansion coefficients 
of within the subspace B\ in powers of A: 



/' 

Jo 



|Aoi^ A ylio|| dto < A n a n (t). (97) 



n=0 



Theorem 7.2 Suppose that 

/•oo 

/ \\ A 01 U t0 A 10 \\ dt <oo. (98) 
Jo 

Suppose that 

a n (t) < c n \t\ n l 2 (99) 

for all t € M and n > 1, where the series X^^Li c n zU has infinite radius of 
convergence. Suppose also that for some e > 0, d n , and all t > 

a n (t) < d n \t\ n / 2 -" (100) 



Then the conditions of Theorem 7.1 are satisfied, namely, there exists some 



< c < oo such that for every r > 

r-\- 2 r 



f'A T 

/ \\AoiU*Aio\\dx < c (101) 
Jo 



uniformly on |A| < 1, and also, for every < r < oo, 



r\- 2 T 

lim / WAq-iOU* -U x )Ava\\dx = Q. (102) 
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,.\ 2 r OO oo 

/ \\A 0l U*A 10 \\ dx<Y, A n a n (A- 2 r) < c n \f\ n l 2 (103) 

n r> 



Proof. By expanding U£ in a A power series, one obtains 

'° n=0 n=0 

which converges for any r. Similarly, 

„A- 2 r oo oo 

/ \\A Ql (U* - U X )A W \\ dx < ^A n a n (A- 2 r) < ^ d n |A| 2e |r|"/ 2 . (104) 

^° n=l n=l 

For |A| < 1 the series is dominated by the convergent Y^n=i c n \f\ n / 2 , and 
each term of the series goes to zero when A — > 0, completing the proof. □ 
Clearly, these conditions refer to the decay properties of the n-point cor- 
relation functions a n (t), and show that the hypotheses ([50]) and ([5T]) are 
satisfied provided information flows fast enough from the subsystem £> to 
the remaining degrees of freedom B\. It is important to note that these 
conditions however only refer to the possibility to perform a semigroup ap- 
proximation, thereby eliminating in markovian fashion the memory kernel, 
and thus only account for irreversibility: they are not directly linked to ob- 
taining a dissipative process, or at least, they are not sufficient. To account 
for dissipation, we need the results of the next two sections. 



8 Dynamical Time Averaging Map 

Our second main result will be to pick, among the possible choices just found 
for the semigroup generator, the most symmetric, and perform a "dynami- 
cal" temporal average, which will always be well defined, no matter which 
are the spectral properties of Zq or the dimensions of Bq. The term "dynam- 
ical" here means that we shall find an alternative averaging map, different 
from the one introduced in [TJ, that depends on the coupling constant, to 
remove the singularities that so severely limit the usual time average. The 
final and remarkably symmetric form of the dynamically averaged opera- 
tor will then be shown to be again compatible with the exact evolution. 
But what is most important, it will be shown to generate a Contraction 
Semigroup in the next section. 

In [TJ a spectral averaging is introduced: for an operator K : Bq — > Bo, 
we put 

K * = 1™ r dq U -* KU * ( 105 ) 

1 ^ + 00 Zl J __rp 



23 



whenever the right hand side is defined. It is important to note here that 
one has to take the limit T — * +00 in order to perform the required time 
averaging (this in fact allows to diagonalize K in the finite dimensional case) . 
Then Davies shows in [2] that if Bq is finite dimensional, then the operation 
\ is well defined and for every r > 



lim < sup 



e (Z a +X 2 K)t _ e (Z +\ 2 K*)t 



(106) 



In the important example studied in [T] , the author finds a completely pos- 
itive dynamics p3] for a finite dimensional system coupled to a heat bath 
using the operator Kjj, where is defined in (|13p . However, in [2] the 
author shows that the temporal average introduce in [1] "was something of a 
red herring" . On the other hand the evolution generated by the unaveraged 
Ko is no more positive in general, so the role of the time average becomes 
instead somewhat important for positivity. Unfortunately, as said before, 
the time average t| is not generally defined when the system hamiltonian 
has continuous spectrum. Since in this work we would like to remain in the 
contest of a generic Banach space, we won't study positivity, but we shall 
focus on more general dissipative properties of the generators. 

Thanks to Theorem 17.11 we are now ready to introduce a new type of 
temporal averaging, that will scale with the coupling constant A and will 
always be well defined (except possibly the singular and uninteresting case 
A = 0). 

We start choosing a = 0, which seems a rather symmetrical case. Now 
let's fix some positive T, and note that 



00 



K {0 , q ,T) = U q \j dxe-^) 2 / T2 U*A iU x A 10 U*\U- q 







U q K 0fiiT U- q . (107) 



The idea is now the following: since in the limit of small coupling A — > 
the system transition time goes to T(A) — > +00, we could use it to diago- 
nalise the generator -K"(o,o,T(A)) by just performing a T(A)-dependent gaus- 
sian integration in the g-variable. This would allow the time average of 
^(o,o,t(a)) to depend on A (hence the name "dynamical"), and to be well 
defined everytime A 7^ 0. To this purpose we give the following 
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Definition 9 For any real positive T > put 

1 f°° f'°° / \2 

= / dqe-i 2 '? 2 / dxe-(V l T U* +q A 01 U x A 10 U^ q 

V^T J-oo Jo 2 2 

1 f°° *f f* 1 '1 

= / dhe-vr* Aoxfa) dt 2 e'^ A w (t 2 ) (108) 

where we have denoted Aij{t) = U-tAijUf. 

Note that in the last line we have changed variable according to 

j * = X/ / " 9 (109) 
[ t 2 = -x/2 -q 

In passing, we observe that Kt so defined can be easily written in more 
physical terms as 

^ r+co r+oo t?+t| 

K T = -_-/ dh dt 2 e~^~ T[A 01 A w ](t 1 ,t 2 ) (110) 

2V 7rJ J-oo J-oo 

where we have introduced the time ordering 

T[AoiAi ](h,t 2 ) = Aoxit^AwihWh-t^+Aotit^AtoihWh-h), (111) 

9 being the Heaviside step function. Indeed, this expression closely resem- 
bles well known von Neumann and Dyson series expansion for the unitary 
evolution operator of a closed quantum system (see for example [17J), and 
to our opinion could offer a great help in understanding how the Nakajima- 
Prigogine-Resibois-Zwanzig master equation ([8]) could be correctly approx- 
imated beyond second ordeiH. 



Theorem 9.1 Suppose that is a one-parameter group of isometries. 



Suppose that there exists some < c < oo such that for every r > 

-A" 2 



M T 

/ \\A 01 U^A w \\dx <c (112) 
J o 



'0 

uniformly on |A| < 1. Suppose also that for every < r < oo 



f'A T 

lim / \\A Ql (U$ -U x )A l0 \\dx = Q . (113) 
J Q 



Another important explicit form for Kt will be addressed in the proof of the next 
section 
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Let T 6 C([— 1,1], R) be a real valued positive continuous function on the 
interval [—1,1], such that 

T(\)~\\\-*f, A~0 (114) 

for some real positive reference time T > and scaling < £ < 2. Denote 
with 

W t x = exp{(Z + AA)0 + A 2 K T(A) )t} (115) 

i/ie associated semigroup on Bo. 
Then for every T > 



lim { sup \\W t x - W t x \\ } = 0. (116) 



0<t<\- 2 T 



Proof. We shall borrow most part of the proof of Theorem (|7.ip . Ac- 
cordingly, we should denote for example with Ti-i^xq) the operator defined in 
Theorem (|7.1j) as W(ao<j)) an d so on. 

Define 

n ij) - - f°°dae' q2/TW2 H U) (117) 

This is obviously Volterra, being an integral of Volterra operators. A closer 
inspection soon reveals that also 

llw?ll -^/*^ /T<A,2ll <» ll - c7 - <118) 

(i) 

as for each j, W^', is bounded by c uniformly on q and 

We proceed on the very same lines of Theorem (|7.ip : the proof that 

lim \\n { ^ - H { i' l) || = (120) 

for j = 1, 2, and j = 6 are in fact identical to that of the foretold Theorem, 
as, for those values for j, one has 

i kn oii w (i)- w ( ( r 9 ) ) ii =o ( m ) 
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uniformly on q. Property (|119p can be exploited to state that 



lim \\H 

A^O 



U) 

A 



n 



(j-i) i 



for j = 3 and j = 5. In fact, for these values of j, we can estimate 



pi 



Ti 



0-1) i 



< 



1 



0FT(A) 



dqe^ T ^' || W 



C?) 



n 



i 



(A, 9 ) ""(A,?) 



1221 



(123) 



Now, the norm in the integrand goes to zero as ~ cf\\) + X 2 q, with 
cf\\) —* uniformly on g as A — > 0, as already noted in the proof of 
Theorem (|7,ip . so the whole integral goes to zero as A — > precisely because 
T(A) scales with |A| - ^ and £ < 2 by our hypothesis. 



It remains to show that if = £ n ( H% ' ) 6, for some initial condition 
b £ Bq, then 



,0) 



lim 



p(4) #(3) 



' A 



A 



. 



(124) 



Proceeding according to Lemma [7J we estimate 



ll/f ~f?\ 



Now we compute 



< EE 



tiu ( n ~ l » 

< 4 >Y , 



-9 2 /T(A) 2 



(3) 



(4) 



(A?) t (A«) 



(125) 



(3) 



(-1) 



(Ag) "(Aflr) 



(i) 



5(A,t,<z) 



H 



("1) 



(cr-A 2 (g-a;/2))- 



(i) 



cr 



1 



V^FT(A) 



< ( K } 



dg ' e - 9 ' 2 /T(A)^ yy ( 

C(A,r,g) 

b 



A 2 ((7-x/2)) - 



< ( K> 



i-i 



\fnT(\) 



( x /2y 



dadxe W K q (X,a,x) 

V(X,T,q) 



ci<Tidxi A A (cTi,xi) e W if^ (A, cti, xi) 



(4) 



/-I 



(ax) (126) 



S(\,q,q' ,a,x) 



27 



where we have put A A = Xi ~ X21 Xi being the characteristic function of 
V(X, a — X 2 (q — x/2) ,q'), x X the characteristic function of V(\,a, q'), and 
we have defined 

S(X, q, q', a, x) = V(X, a - X 2 (q - x/2) , q') U V(X, a, q') . (127) 
Passing to the norms we obtain 

x sup [[ dadxe~ {x/2)2/TW2 k(x) E qql (X,a,x) (128) 

0<T<T J J 

V(\,T,q) 

where the slight modification of the related definition of H g (A, a, x) in (jUJ) 
is given by 

E qq i(X, a, x) = Jj do\dx\ \A x (a\, xi)| k{x\) . (129) 

S(X,q,q' ,<r,x) 

Proceeding on the same lines as in Lemma[7]we find the asymptotic behavior 
sup / / dadx e-^^ 2 / T(x ^k{x) E qq ,{X, a, x) 

0<T<T J J 

V(X,r,q) 

POO 

<t sup / dxe- {x/2)2,TW2 k(x) E qq ,(X,a,x) 

0<a<-Jo 

<T{c 1 (X) + c 2 X 2 \q\+c 3 X 2 \q'\} (130) 

where ci(A) —* uniformly on q and q' . This, plus the fact that £ < 2 by 
our hypothesis, allows us conclude that Eq. (|124p holds, as can be seen by 
putting result ()130j) into (|128j) . and then back into (|125p . and by using the 
dominated convergence theorem. 

Collecting the results as in Theorem (|7.ip concludes the proof. □ 

To make contact with the definition of the time average proposed in [1] , 
we give the following 

Proposition 10 Let Bq be finite dimensional, and Aqq = 0. Then for b £ 
Bo and every r > 

lim{ sup \\e( z °+ x2K TwH-e( Zo+x2K >b\\\ =0 (131) 

0<i<A" 2 r 
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Proof. Let 

Zq = iUgQa (132) 
a 

be the spectral decomposition of Zq, with all u; a 's distinct and real, and 
compute 

lim K T = lim — / dqe-^/ T2 y2e l ^-^Q a K m+oo) Q 

= ^ QaKoQa 

a 

T^+oo 2T J_ T ^ 
= Urn / dqU^ q K D U q (133) 

T— >+oo Zi 

which clearly shows that the time average in pQ coincides with our dynamical 
one, in the weak-coupling limit A — ► 0, that is, recalling that lim^oT^) = 

+oo, 

K) D = lmiK T{X) (134) 

and the averaging map \ is well defined because Bo is finite dimensional. 
Then a slight and straightforward modification of Theorem 1.4 in [2 J proves 
the stated result. □ 



• We stress that, contrary to Kjj, our time average K T rx\ is always well 
defined, irrespective of Zq spectral properties, for any nonzero value 
of the coupling constant A. 

• The statement of the proposition would remain unchanged if we relax 
the hypothesis that Bq is finite dimensional, and only assume that Zq 
has discrete spectrum. 

10.1 Comments To The Theorem 

The last line in Definition [9] indicates a high level of symmetry and simplicity, 
so that one could ask whether such a form would be "unique" to some extent. 
With respect to this, we remark here that there is nothing peculiar in the 
choice a = that led us to the dynamical time average in Definition [9l apart 
from simplicity in the definitions and in the proofs involved. In fact, one 
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could equally well proceed along the following lines. In equation (jSJ) change 
variable according to (|16|) as in Theorem 17.11 Proceed along the lines of 
Theorem 17.11 but use 

exp{-h(x, q) 2 + t 2 (x, qf - 2q 2 /2T(X) 2 } (135) 

instead of exp{ — (x/2) 2 /T(A) 2 } as gaussian smoothing for the kernel in ([75]) , 
where 

h(x,q) = (a + l/2)x-q , . 

t 2 (x,q) = (a-l/2)x-q ' 

There is no difficulty to proceed as in Theorem 17. II to show that (under the 
same hypotheses) the semigroup = exp{(Zo + A^oo + )?K(os,q,T(\)))i} 
satisfies 

lim \ sup ||W t A - W t X \\ \ = 0, (137) 

[o<i<A- 2 r J 

for every r > 0, with the modified version 



K (a, g ,T(x)) = dxe U_ {a+ r )x+q A m U x A w U {a _, )x _ q 



oc 



dxe 4u(*i0c,g))>M*2 0c,?)) ■ (138) 

Then as in Theorem 19.11 one could perform a gaussian integration on the q 
variable to obtain 



1 /"CO 



l r 00 *i r* 1 «1 

= / dtie ^WAoi(ti) / dt 2 e ^4(#9) 

exactly as in Definition [9j This shows that the choice a = is merely dic- 
tated by the simplicity of the proofs involved, and furnishes an argument for 
the uniqueness of expression [9j One could argue that the choice of a gaussian 
to smooth the involved kernels is somewhat arbitrary, but, as commented in 
|13j . the gaussian is the only distribution $ with the factorization property 



$(g)$(s/2) = $(ii/V2)$(t 2 /V2) (140) 

for t\ = q + x/2 and t 2 = q — x/2. As the remaining degree of arbitrariness 
is concerned, the transition time T(A), we have already observed that it is 
actually fixed by the obvious and natural choice T(A) = 1 1 1 1 1 . 
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11 A Contraction Semigroup 

To establish the extent to which is unique, we recall that % is a con- 
traction semigroup if \\%\\ < 1 for all t > 0. 

Before stating our main result of this section, we report here for com- 
pleteness the Hille-Yosida Theorem [18) . an important part of which will be 
used throughout in our statements (although the operators involved will be 
bounded). We denote with R(S) the range of an operator S. 

Theorem 11.1 (Hille-Yosida theorem) Let S be an operator on the Ba- 
nach space B. For T = £>* or T = B* the following conditions are equivalent: 

1. S is the infinitesimal generator of a a (B, J-) -continuous semigroup of 
contractions T ; 

2. S is a {B^T)- densely defined, and a(B,J r )-a(B, J-")-closed. For a > 

||(l-a5)6||>||6||, b€D(S), (141) 
and for some a > 

R(l -aS) = B (142) 

If these conditions are satisfied, then the semigroup defined in terms of S by 
either of the limits 

%b = limexp^l-eS)- 1 ^ 

= lim (1 - St/n) n b (143) 

n— >oo 

where the exponential of the bounded tS(l — e5) _1 is defined by power se- 
ries expansion. The limits exist in the o~(B,J-) topology, uniformly for t in 
compacts, and in norm if b £ D(S) is in the norm closure of D(S). 

We shall also need the following: 

Lemma 12 Let A and B be bounded generators of one-parameter groups of 
isometries on a Banach space B. 

Then C = [A, B] is the generator of a one-parameter group of isometries 
on B. 

Proof. Define 

F t = exp{^t} exp-t^t} exp{-^t} exp{-5t} t G R (144) 
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Clearly Ft is an isometry for every t E R. Now for t € R \ {0} define 

C n (t) = n 2 t~\F t/n - 1). (145) 

Then C n {t) generates a one-parameter group of isometries. To show this we 
take a > 0, b £ B, and prove inequality (|14ip by computing 

||(1 - a(F t/n - l))b\\ = (1 + a)||(l - -^F t/n )6|| 
' 1 + a ' 



> (1 + a) 



" "F t /nfe|| 



1 + a 



> IN, (146) 



as H&ll > a/(l + a) because F t / n is an isometry and a/(l + a) < 

1. It can readily be seen that all the other hypotheses of the Hille-Yosida 
Theorem are satisfied (recall that A and B are supposed to be bounded). But 
proceding in the same way, inequality (I14ip can be seen to hold also for a < 
0, so that C n (T) actually generates a one-parameter group of isometries [19]. 

Now a simple calculation [21] shows that lim n ^ 00 C n (t)b = Cb for ev- 
ery b G B. But the inequality (|14ip . together with the validity of all the 
other conditions of the Hille-Yosida Theorem, passes to the limit, due to the 
boundedness of the involved operators, thus proving the Lemma. □ 
We can now state the main result of this section: 

Theorem 12.1 If \\P \\ = 1, then W t x is a contraction semigroup on Bq, 
for all real X. 

Proof. First of all X x is a one-parameter group of isometries because of 
Lemma [3l Now because of the Trotter product formula [20], one has 

P W t X = p oe (Zo+XA OQ +\iK T(x) )t = Hm ( X X e ^K T(x) t/ny (147) 

so that 

\\PoWf\\ < sup HPoe^W^ir, (148) 

n 

and the theorem would follow if Kt would generate a contraction semigroup 
on Bq, for all T > 0. To show this is indeed the case, we name 



m = J^e-^U. t AU t (149) 



and denote as usual &ij(t) = Pi&(t)Pj. Then from 

-oo p—oc r+oo 



/+oo r+oo r+oo rti 

dh / dt 2 $ l(*l)$10(*2) " / dt X / rft 2 *0l(*2)*10(tl) 
-oo J —oo J —oo J —oo 

/+oo nt\ 
dh / dt 2 $oi(*i)$io(*2) (150) 
-oo J —oo 
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we sum term by term to obtain 

R t=~ d*i$oi(*i)/ dt 2 $io(t 2 )+- dh dt 2 [$ i(*i),$io(*2)] 

(151) 

By naming C = J*_ °° (it <J>(t) and Coo = ^ 3>oo(*)> we introduce 

= ^ (C " Coo) 2 + g / *i / *2 [(<*> - *oo)(*i), (* - *oo)(*2)] (152) 

where (<E> - 3>oo)(*) = - $oo(*)- It follows easily that 

K T = P K T P . (153) 

Now again through the Trotter product formula, the theorem would follow 
from the fact that Kt generates a contraction semigroup on B. Indeed, 

P e KTt = lim P {1 + P K T P t/n + 0(n" 2 )} n 

= lim {P e RTt/n P } n , (154) 
so that, since ||Pq|| = 1, we obtain 

\\Poe Krt \\ < \\e KTt \\ (155) 

for all t > 0. So we shall prove that Kt generates a contraction semigroup 
on B by showing that each of the two terms in (1152|) does. 

In order to do that, we note that both C and Coo (and thus also C — Coo) 
generate a one-parameter group of isometries on B, as can be easily seen 
through repeated use of the Trotter formula for both positive and negative 
t, recalling we assume A and ^4oo generate isometries, and Ut are isometries. 
In fact, recalling C is bounded (as Coo), and that $(t) — ► for t — > ±oo, 
with obvious notation one has 



e ct = lim hm{TTe A ^ fe)t/n } n (156) 



and for all subdivisions {Aj} and natural n 

||{JJ e A,-*fe)t/n|n|| < "Q UgA^^Jt/niin < ^ ( 157 ) 

j 3 

as e*(* is a contraction on B for all real t' and t. This in turn follows 
because &(t') satisfies all the hypotheses of the Hille-Yosida theorem pjj] as, 
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being bounded for every finite T > 0, it is in particular weakly closed and 
densely denned, and for all real a > \\A\\ > ||<3?(i')|| the range R((31 — &{t')) 
of (31 — &(t') is the entire space B (see propositions 3.1.1 and 3.1.6 in [T9]). 
Moreover, for all real a and b £ B we have 

||(1 - a*V))H > 11(1 - -=A)U t b\\ > \\b\\ (158) 

because A generates isometries and U±t are isometries. Then the Hille- 
Yosida Theorem asserts that exp{<3? (£')£} is a one-parameter contraction 
semigroup (in the t variable) for every real t' . Since the inequalities above 
hold for all real a, Corollary 3.1.19 of [19] further implies that exp{&(t')t} 
(and hence also exp{Ci}) is a one-parameter group of isometries on B. The 
same arguments imply that Coo generates a one-parameter group of isome- 
tries on B. 

Now the first term in (I152p generates a weakly-continuous (and hence 
strongly continuous) semigroup of contractions through the Hille-Yosida 
Theorem [19] as, being bounded for every finite T > 0, it is weakly closed 
and densely defined, and for all real a > ||-Kr|| the range R(al — (C — Coo) 2 ) 
of al — (C — Coo) 2 is the entire space B (see proposition 3.1.6 in [19]). More- 
over, by writing 1 - a(C - C 00 ) 2 = (1 - V^{C - C 00 ))(l + V^{C - Coo)) 
for a > 0, and using the fact that C — Coo generates isometries, we see that 

||(1 _ a{C - C m ?)b\\ > ||(1 + v^(C - C 00 ))6|| > H&H (159) 

for any a > and b £ B. 

To treat the remaining term in (|152p . we denote Dt = ($ — ^oo)^)) s ° 
that by the Hille-Yosida Theorem D t * generates a one-parameter group of 
isometries for any real t'. Now, according to Lemma [T2| [Dt t , Dt 2 ] generates 
a one-parameter group of isometries for every t\ and t<i- The integrals on 
such t\ and t<i can again easily be shown, as before, to generate isometries 
by the Trotter formula (they generate contractions for all t 6 R). 

This completes the proof. □ 



12.1 Comments to the Theorem 

• It is evident that none of Ki a ^j>) generates contractions, so this the- 
orem actually establishes a "uniqueness result" in the class of all the 
possible semigroup approximation of the exact projected dynamics (jSJ). 
Moreover, having a contraction semigroup at disposal is very impor- 
tant for physical applications, where one is willing to study the limit 
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dynamics at all times, including steady states. It is well known indeed 
that a Quantum Dynamical Semigroup in a von Neumann algebra [2] 
is, in particular, a contraction semigroup. 

• We have shown that our dynamical time coarse-graining procedure 
does lead to a dissipative generator for the irreversible dynamics. Phys- 
ically one could say that a certain scale is fixed (the transition time), 
in which the subsystem is not able to distinguish the details of the evo- 
lution. Performing a time average on that scale allows one to eliminate 
the exponentially growing "off-diagonal terms" in the generator, thus 
obtaining a genuine dissipative process. Probably this could be ex- 
plained by saying that K^ a q T ^ "sees" pairs of transitions A m and A w 
instantaneously, at two specific and different instants inside the tem- 
poral transition window/scale, resulting in a subsystem internal transi- 
tion which is strongly asymmetric (and singular) in time. Through the 
dynamical time average Kj>(\) instead, subsystem internal transitions 
take place smoothly and homogeneously in time. 

• Given that Kt = PqKtPq with Kt in (|152p is the sum of two terms, 
the first term is dissipative, whereas the other one is a second order 
energy renormalisation, generating a one-parameter group of isome- 
tries. 

13 Summary and Conclusion 

We have considered the class of quantum mechanical master equations for 
Physical Subsystems in the weak coupling limit. In the first part we have 
taken this to be a perturbed one-parameter group of isometries projected 
on a Banach subspace. In that limit, we have shown that the memory terms 
of the integral equation for the evolution operator can be approximated in 
a variety of different ways: the resulting class of one-parameter semigroups 
includes the previous literature. Then through the introduction of a dynam- 
ical time averaging map we have found explicitly a generator that i) is able 
to correctly approximate the exact projected dynamics in the weak-coupling 
limit, ii) is always well defined, irrespective of the subsystem spectral prop- 
erties and dimensions, Hi) accounts also for first order contributions, iv 
boils down to previous literature results in case of discrete spectrum, and 
more importantly v) gives rise to a contraction semigroup on the projected 
Banach subspace. 

Our results are of very general nature, and open the way to the study of 
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a variety of different extensions (beyond second order, time dependent free 
dynamics, to name a few) and physical applications in the field of operator 
algebras, also greatly linked to nowadays technologies. 
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